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10. Prediction of Large Negative 
Shaded-Side Spacecraft Potentials 


1 M. L PrQkdftenko dnd }. G. Ldfrdmbois^ 
Physics Dsporfnteftt 
Y0ik University 
T0rbittb« Canada 


Ai>stract 


A calculation by Knott, for the floating potential of a spherically symmetric 
synchronous -altitude satellite in eclipse, has been adapted to pro\tlde simple cal- 
culations of upper bounds on negative potentials vrhiCH may be achieved by elec- 
trically isolated shaded surfaces on spacecraft in sunlight. To investigate geomet- 
rical effects, we have replaced Knott's use of the orbit -limited ion current expres- 
sion for a sphere, by that for an infinite cylinder. Large (‘^60 percent) increases 
in predicted negative shaded-side potentials are obtained as a result. To investi- 
gate "effecUve -potential barrier" or "angular -momentum selection" effects due to 
the presence of leSS-negative sUnlit-Side Or adjacent -Surface potentials, we have 
also replaced these expressions by the ion random current. Which is a lower bound 
for convex surfaces when such effects become very severe, f^urther large in- 
creases in predicted negative potentials are obtained, amounting to a doubling in 
some cases. Depending on Surface properties and incident energy distributions, 
values exceeding -20 kV are now predicted, in good agreement with ATS-6 obser- 
vations of potenuais reaching -19 kV, aS reported by Whipple. For isolated Sur- 
faces in shaded cavities, even more negative values may be reached, m some 
conditions, two distinct, floating potentials exist, leading to the possibility t 
"bifurcation phenomena" in real Situations. 
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The performance of many satdliites in synchrortous orbit has been degraded 
by aitomalous evertts which include frettueitt spurious spacecraft commands and in 
some cases permanent damage. These events irtvariably appear to involve elec- 
trical discharges caused by differential charging of spacecraft surfaces to large 
relative potentials. The latter condition in turn is known to result from the rela- 
tively large average energies (Up to a few' keV) of the charged particle environment 
at synchronous orbit altitude, particularly in disturbed magnetoSpheric conditions. 
Since photoeiectron emission from sunlit surfaces tends to compensate for incident 
electron fluxed, thereby holding sunlit -side surface potentials close to space 
potential in rnost cases, an estimate of differential charging magnitudes can be 
obtained by simply calculating floating potentials of electrically isolated shaded 
surfaces, relative to space potential. In the present work, we have at "^mpted to 
obtain Upper bounds on such potentials, which in cases of interest are usually 
highly negative, because these bounds constitute "worst cases" for design purposes, 
and also because unlike more exact calculations, they can be obtained from simple 
current balance calculations. Furthermore, it is sufficient to consider local cur- 
rent balance only, because this corresponds to an electrically isolated surface 
element, w'hich is also a "worst case” for differential charging. To calculate these 
bounds, we have extended a calculation by Knott, ^ of the floating potehtial of a 
spherically symmetric svTichronoUs -altitude satellite in eclipse. To investigate 
geometrical effects, we have replaced Knott's use of the Mott -Smith and Langmuir^ 
orbit -limited curCent expression for collection of Maxwell i an ions by a unipotential 
sphere, by the corresponding expression for an infinite cylinder; both expressions 
have been shewn to be upper bounds for collisionless ion collection as a function 
of lotal surface potehtial, for three'* arid two-dimensional collectors, respectively, 
regardless of collector shape, sheath potential, or potential of other parts of the 
collector. This replacement causes a large decrease in ioti collection and a 
correspondingly large increase in negative shaded-side floating potentials (Sec- 
tion 3). Another important ion-current restriction mav be caused bv '’effective- 
potential barrier’"'^' • ^ or ’'angular-momentum selection” effects, in which the 
presence of less-negative sUnlit-side potentials produces dipole and higher mom- 
ents in the sheath potential, ^ causing steepening and contraction of the potential 
well surrounding the shadad side (Figure 1). A similar steepening effect will also 
occur if an isolated shaded surface element is surrounded by adjacent shaded 
surfaces which for any reason have less negative potentials (Figure 2). The most 
extreme possibility would be a potential profile which was equal to space potential 
almost to the spacecraft surface, then fell discontinuously to surface potential. 

This limit would correspond to a "planar sheath” situation in which the ion collection 



Figure 1. St^ept^niiig of Shadfed 
Side Potential Profile fof a 
Spacecraft with an Insulated Sur- 
face, after Fahleson” 


Figure 2. Conductive Spacecraft with 
Shaded Isolated Surface Patch 


on any shaded convex surface would be given by just the ion random flux. This 
amounts to a further ion cUrreht restriction which produces even larger increases 
in negative shaded-^ide floating potentials (Section 3). This situation corresponds 
to a velocity -space cutoff boundary for incident ions which is “one-dimensional;" 
the cutoff boundaries corresponding to spherical and infinite cylindrical collectors 
respect lively, “tht'ee-dimeftsiortai" and **two-dihienslonar’^ (Section 2). 

We also show (Section 3) that if shaded cavities containing isolated surfaces 
exist on a spacecraft, negative potentials ort such surfaces may surpass even these 
predictions. In some Cases, more than one possible floating potential results 
from the calculation; this implies the possibility of ’’bifurcation phenomena” in 
\Vhlch adjacent isolated surfaces of the same material may arrive at different 
floating potentials as a result of differences in their charging histories (Section 3). 

We have also modified Knott’s calculation in another way, by including currents 
due to electron backs battering (Section 2). These currents ^ill tend to decrease 
net electron collection, thereby making floating potentials less negative than other- 
u^lse (Section 3). A process not included by either Knott or ourselves is secondary 
electron emission due to ion impacts; this will also tend to make floating potentials 
less negative. 
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2. TllWtllV 

The ambient eleetfon energy nutribntlena uaed In the |,r. aent «ork are a 
model dnlet-tlme apeetmtn (Knelt.' Figure t) and a modei disturbed spocirun. 
(Kftott Figure 2b) ba.sed on meajiUrementaS by Shield and ^Yank, and DeFoie&t an 
McIlV^dn, ^ respectively. Both of these dlstrlbutloas. and al.so the ambient ion 
distribution, are assumed isotropic. The disturbed spectrum was 
three used by KhOtt because it has a higher- average electron ertfiCgy ('10 keV) 
than the others, in using it. we have chahged it as follows: in the --gy rang- 

0 5 keV ^ E 10 keV and 10 keV E s 40 keV. we have replaced Knott s differen 
» ..V. Knr h y 10® and \2 X 10® E”^^^ eleCtrons/cm sec sr 

tlal energy spectrum by n'2 X 10 E anoN.«Aiv u, 

keV. respectively, where E Is energy. These r-’ .ionsareslmplerthantho.se 
indicated by Knott, and they also bring the model spectrum into closer 
with the data on which it is based. We therefore believe that they may have been 
the ones actually used by Knott, and that the corresponding parts of Mgure 2b m 
his paper may be Incorrectly plotted. For any spacecraft surface having a nega- 
tive potential d . < 0. or for a three-dimensional (for example, spherical) surface 
having d. > ‘he orbit -limited BuX (particle current density) of ambient elec- 
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trons is given by: 


•I. ■ / "" 

max(0. -ed^,) 


edg/E)(djp^/dE) dE 


( 1 ) 


where e is magnitude of unit electronic charge, d^ is local surface potential 
dJ. /dE is the ambient energy-differential flux incident on one side of an arbitrar- 
IW oriented surface element, and v„ is the velocity component normal to the same 
surface element. dJ^^/dE is ir times the energy-differential flux per 
used by Knott. ^ and Is given m terms of the ambient electron velocity distribution 
f - d® PJ /d\ by the relation dJ^^/dE = 2»rf£/m2. where m^ is electron mass and 
N is ambient ion or electron number density. Since f is isotropic, f s f(E). The 
factor (1 + ed /E) in Eq. (1) appears to have been neglected by Knott, and may 
account for some minor discrepancies between his results and ours (Section 3). 

The presence of this factor leads to a divergent integration in Eq. (1) 'f dg > 0. 
unless dJ.^/dE - 0 as E 0. that is. f(fi) remains finite as E - 0. This 
that the dlrterentlal fluxes in Kriott<s Figures 1-3 must approach zero linearly With 
E at E values smaller than those shown IH these figures. In the present work we 
have instead used a less-realistic sharp cutoff at i eV; this may slightly affect our 
results for positive noatirtg potentials in Section 3. We have also Introduced a 
sharp upper cutoff at 50 keV in (he quiet -time spectrum, also in order to avoid a 
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dlver^ortt intejjratldn when ralculatlng avei’ag^ enei^gy for use In backscatterlng 
ealctilations (see below). 

In order to obtain the orblt-llmited electron flux expeossion for an arbitrary 
cyllndrleal collector, tho lower integration lltnlt in fig. (l) must be replaced by 

the two-dimensional velocity. space cutoff boundary IS max (0, - eo ) whore 

h is the total energy of transverse motioh ^ n^^^ . - c«. and we have 

chosen a .. coondinale perpendicular to the cylinder cro.^s-soetion. if o - 0 this 
complicates the integration in Eq. (l). which may then be done in eithcr'of two 
ways. The first | Eaframboise and Parker. Eq. (,i)| is to convert Eq. ( 1 ) into a, 
in egration using cylindrical coordinates in velocity space. This method has the 
disadvantage that the velocity distribution must be integrated over v in order to 
convert it into a distribution of transverse velocities. Art alternate method=^- » is 
as followd we Choose rectangular coordinates (v^. v v^) in velocity space, such 
that v„ ,s the velocity component in the inward normal direction at the collector 
surface. The x and v^ become tangential coordinates, with v. in the plane of the 
cylinder cross-section. We Ihen transform to spherical coordinates (v. e v) with 
as polar axis. Then: v cos v v sin o cos and v, v sin 0 sin 

FOb > 0. Eq. ( 1 ) is then replaced by: t » t'- 




E-^ 0^7t/2 0~7t/2 

f. / f(£)(v siri 0 cos 0) 

E 0 «; -ff/2 0 ArcsirtIe0^/(E + eOjj)l 

2 

X (v sin 0 dv dfl d(/,) 


h\ 


Arc sin 


E + ^ “FT 


^ J V E j 


dj 

dE . (2) 


extra 0 ^ 1 " a ' Eq. (2) contains an 

extra energy-dependent weighting factor, which arises from integration of v 

over the fractional solid angle over which ambient eleetbons can beach the coHec 
tor at eabh energy, 

advantageous in ' tatnlng the one-dimensional 
(Sec ion 1) orblt-llmlted Hux expression. IH this case, the lower limit in Eq. (1) 
must be replaced by: E„ max(0. - ee,). where E„ » m„ vf - ee. This til 


we ti*anBform <v^, Vj, to spherieal doordiftatea (v, 6 , «/) with aa polaf akia. 
for > d, Eti» (1) ts now reiil&oed bys 

Epo d=Afcc*oaledg/(E+eig)) 

Jg “ J J f(E)(v Ooa d)(2ffv^ silt d dv dd) 

E*0 d=0 



<3) 


independently of collector potential, as expected. 

The corresponding expressions for ion flbx are simpler because the ions— 
are assumed to be Maxwellian. Corresponding to the three-, two-, and one- 
dimensional velocity-space cutoffs described above, we obtain, respectively, ^ 
for~ion-attractlng surface potentials Xjg > Ot 




= J 


to 


<1 ^ 

(2(Xjg/ff)^/^ + expix^g) *rfc 

(d) 


(4) 

(5) 

( 6 ) 


Where x^g = -edg/kT^, k is Boltzmann's constant, Tj is tort temperature and J. is 
the ion random flux N^ikT-./2ffm.)^^^. J'or lon-retardirtg surface ptotentials 
Xjg < 0, we obtain; 


Ji = exp (xig) . (7) 

We have assumed^ -that Tj = 1 keV, and that the random ion to electron flux 
ratio = i).025. 

Pot the secondary electron fractir nal yield 6(E), we have used, following 
Knott, ^ the relation of SterngUss;^^ 

= . («) 

We have used the same selection of surface materials (Section 3) as that 
appvsaring in Table 1 of Knott, ^ for Which the values of and E ^ were 
obtained froih Glbbdtis ^ and ti&ch^nberg and BraUer* 

The process of electron backscatterlrtg, which was not included in Knotts 
calculations, becomes important at incident electron klhetlc energies larger than 


those for which secondary emission Is dominant. For th^ backscattored . lectron 

fractional yield r,. ^ve havd fitted the results of Stcrnfilass^"^ with a I'Mati.ih of the 
forfn: 

n(E) AE + t C 

where the coefficients A, B. and C are functions of the atomic number Z of the 
surface material. WC-havC evaluated A, B. and C for each surface material con- 
sidered (Section 3) by substituting StCrnglass' values of n at l, 3. and 5 keV, into 
Eq. (9). m all cases, n is a very slowly varying function of E. For compound 
surface materials, we have assumed that each atomic constituent contributes an 
independent backscattered nux proportional to its relative concentration. There 
exist more recent measurements ot which give generally larger values than 

tho.-e of Sternglass, especially for electrons having near-tangential incidence. 
However, we have found these results to be too fragmentary for our purposes, and 
we have therefore used Sternglass' results throughout. Presumably we have 
therefore underestimated ij. and our predicted floating potentials in Section 3 will 
therefore be slightly more negative than more realistic corresponding values. 

tthen > 0. not all secondary and backscattered electrons will escape. To 
calculate flux escaping, we assume^^* for ease of calcv’ation. that both secon- 
dary and backscattered electrons are emitted with Maxwellian velocity distributions 
having thermal enfefgies 

^sec " *^^sec ^ ^scat ' *‘^scat + 2 X lO^^Z)!" eV 

Respectively, regardless of the form of the incident velocity distribution. Here. 

E IS the average incident electron energy. We further assume that escape of 
emiUed electrons is orbit-limited, that is. that no barriers of effective poten- 
tial • • or negative barriers of electric potential exist on the shaded side. 
Fahleson has pointed out that such barriers are likely to exist on the suhlit side 
independently of any space-charge effects, if substantial shaded-sunlit differences 
exist in The expressions for the escaping secondary and backscattered fluxes 
'^scc *^ 3 cat ^^>®*’®fore aret 
oe 

•^sec ” y* e<>^/E){cUpy/dE) dE dO) 
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( 11 ) 


•’scat I 

4 ^, 

If 0 <0, Ifv >0, the three-, two-, and one-dimensional cases must be fconsid- 

ei^ed separately. We define e<i> 3 /kT^gg and ^ 

w.e consider only the secondary fluxes; the corresponding results for backscattered 
fluxes may be obtained by replacing 6 by n and hy Xj^^gj throughout. If «ljj >s 
the emitted flux of secondaries, then their velocity distribution at the surface is 
’s = <*’s/^*’^'^e^*'’’’8ec’^ exp ( - 7 m^v^/kT^ ). In the three-dimensional case, 
the cutoff condition for their escape is E - J m^v^ - ee^ > 0. We redefine as 
velocity component in the outward normal direction, and we use spherical coord- 
inates as defined in connection with Eq. (3). We obtain, for the escaping secondary 
flux^ — 


sec 


/ 


''n "’l 



hz L 

E=oo ® 2 

/ / exp (- j e)(2irv^ sih 0 dv dg) 

E-().. 0=0 


= ^sec’ ®*P ^■’^sec’ J ^ e^g/E)(dJg^j/dE) dE 


( 12 ) 


The factor. (1 + X is noteworthy, because it Is specific to three: -dimensional, 
as opposed to planar, sheath geometry. In the two-dimensional case, the cutoff 
condition for excape Is m^lv^ + v^ ) - edg > 0, and the convolution integral for 
Jg contains the extra weighting factor which appears In Eq. (2). It Is convenient 
to use spherical coordinates as defined in connection with Eq. (2). We obtain: 


I 


I 


'tide 




1/2 


exp (X 


1/2, 

!^ec' 


oxp ^Xrtoc' 



K + \ 


-^j 6(K ^ e<!,^) 


V 



dE 


(13) 


12 A ^ 

In the one-dltnenisional case, the escape condition is ^ ^ 

we again use spherical coordinates as defined in connection with E<lt (3), We 
obtain: 


go 

^sec ®’‘P ^■’^scc’ ^ e0^)(dJgy/dE) dE 


<14) 


The floating poterttial(a) nf an isolated shaded surface element Is (are) nov? 
given by the zero(s) of the function: 


"^net " "^e ^ "^sec ^ "^scat 


(15) 


3. RESL'LTS AM) DISr.liSSION 


Table 1 shows floating potential values obtained by numerical solution of the 
ecluation = 0, where is given by Eq. (15) and wb have assumed ^ 0 

irt order to duplicate the physical situation of Kriott, ^ whose results are shown in 
parentheses. We sbe that our results show qualitative but not quantitative agree- 
ment with his. Possible reasons for the disagreement are: (1) Knott appears to 
have solved the current balance equation graphically rather than numerically; 

(2) wherever his solution indicates a floating potential more negative than -3000 V, 
he has listed the result simply as "< -3000 V"; (3) wherever he has obtained a 
positive floating potential, he has listed it simply as " + 5 V" whereas we have 
calculated it using the assumptions made in Section 2 ; (4) as mentioned following 
Eq. (1), his expression for incident electron flux may contain an error. The most 
Important feature of Table 1 is the very large floating potentials which are evident 
in disturbed conditions in the presence of the two- and one -dimensional velocity- 
space cutoffs. The dramatic differences which exist among these results are 
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evidence thal spacecraft geometry and sheath potential shape are important influ = 
oneos in determining floating potentials. It is important to note that as floating 
potential becomes more negative, it also becomes more sensitive to the presence 
of small amounts of high-ener^ electrons. This moans that If a spacecraft should 
encounter conditions that are "more disturbed" than those given by Knott's spec- 
trum 2b, the values in Table I most likely to bo significantly exceeded are those for 
the one-dimensional cutoff. Thir implies that for design purpe les in which worst- 
ease information is desired, it is important to do calculations with tb ■ "most 
disturbed" electron spectra available. 

In obtaining these results, We have made no attempt to calculate the time 
needed to approach the steady-state conditions which they represent. In general, 
the most negative potentials correspond to a balance between the smallest currents, 
and will therefore involve the longest charging times. 

Also evident in Table 1 are situations in which the current -voltage character- 
istic of the surface has three roots. For these to occur, it is necessary that 
^max substantially greater than one, and that the incident spectrum contain a 
sufficient proportion of electrons in the energy range where secondary emission 
is a maximum. The center root never represents a possible floating potential, 
because it is "unstable" in the sense that a small change in surface potential would 
cause a net current collection of a sign which would drive the surface potential 
away from this root to one on either side. A further consequence of such a situa- 
tion is discussed beloW. 

Table 2 includes the further addition of backscattered electron flux (Section 2) 
and therefore represents a more realistic physical situation. In most cases, the 
net effect of backscatter is a moderate reduction of negative floating potentials. 

In some cases, the reduction is large, as in the case of a gold surface exposed to 
the "quiet" spectrum. In several other ca;. es, all associated with the quiet spec- 
trum, backscaitering changes a multiple -roc* to a s*ngle-root situation. As indi- 
cated in Section 2, we have probably underestimated backscattered fluxes, and we 
have also (Section 1) ignored secondary electron emission caused by ion impacts. 
Both of these effects would tend to further reduce negative potentials. However, 
such changes are likely to be small. The results in Table 2 should probably be 
regarded as consistent with observations of po..,'ntials reaching - 19 kV on the 
ATS-G spacecraft, as reported by Whipple. 
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Plgiilfes 3-t sho'w &ui*f6ni -Voltage ehafaeteriatlcs for some bf tRe situations 
In Tables l aAd 2. Figure 3 shows a "typleal" sthgle-robt situatlbn in which 
secondary and backscatter Contributions do not change the general shape of the net 
current Curve. Figure 4 shows the above-mentioned case of gold exposed to the 
quiet spectrutn, In Which the backscatter contribution changes a large predicted 
negative floating potential to a much smaller value. Figure 5 shows a triple-root 
situation. Figure 6 shows the disappearance of a triple-root situation because of 
backscatter. in Figure 7, secondary electron current is sufficient by itself to 
prohibit a negative floating potential. 

We now examine situations which may arise tn the case of spacecraft which 
have shaded cavities containing electrically isolated interior surfaces.^ Figure 8 
shows an idealization of such a spacecraft. We wish to show that the effects of 
surface concavity may Cause ion Collection to be reduced more than net electron-—- 
collection at an interior point sUch as B. relative to an exterior point A; such a 
situation wculd result in floating potentials more negative than those of Table 2. 

To demonstrate this possibility, we first note that in the presence of an isotropic 

q 

ambient plasmk^ inbid^nt fluxes to any surfgice depend only on the locations, in 
velocity space, of the cutoff botihaariC^ IhSide of which the orbits of ambieht par- 
ticles can conhect "from iftfinity*' to the surface. Figiire 8 Shov.^s a- set of the 
associated "cutoff orbits." We see from Figure 8 that the indluded ahgle between 
Cutoff orbits has beCh reduced ih gcing fhorn A to B for ions but hot for electrons, 
for vi^hich orbits tahgehtial tc the surface are shd^h as reaching both A and B. 
Accordingly, the incident ion current contribution for the energy shown will also 
be reduced, but the electron contribution will not. This picture is invalid for 
higher-energy electrons at B, whose orbits are straighter and will have a greater 
tendency to connect back to the interior surfaces of the cavity. Bveh though such 
higher-energy orbits will generally have lower populations than loWer-energy 
orbits, it is not clear whether the relative current reduction at B will be greater 
for ions or for electrons. However, this argument is intended to demonstrate only 
the possibility that the bounds In Table 2 will be e^tceeded. Oh the other hand, this 
possibility will be enhanced by the effects of secondary and backscattered electrons, . 
which will tend to be recollected inside any cavity, rather than escaping Into Space, 
thus tending to increase net electron collection and driving floating potentials more 
negative. This effect will be strongest for backscattered electrons because their 
higher emission energies will cause them to have straighter orbits. To draw firm 
Conclusions will require detailed numerical simulation. An additional feature of 
cavities is thctr generally higher oulgasslng pressures, which will increase any 
tendencies for arcing to occur. More negative floating potentials may also result 
if the ambient elebtroti distribution dontaind beam-like constituents which happeu 
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Flgtire 3, Currant-Voltage Gharacterlatlc for Aliuiilnum in "Quiet" 
Conditions, with a Ohe -Dimensional Velocity-Spade Cutoff. In Figures 
3-7, the zeros of the characteristics are ludicsdted by arrcws 
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Flimre t. Curreht -Voltage Chairacterlstlc for Beryllium -Copper in 
"Quiet" Conditions, with a One -Dimensional Velocity-Space Cutoff 
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SPACECRAFT 


SUNLIGHT 


tCN ORBITS 
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Fiaure B. Spacecraft with Shaded Isolated Cavity, Showing incident Ion 
and Electron Orbits with Energies Close to the Lowest for wMch Collec 
tlbh of Each Species is Possible. The orbits shown are cutoff orbits, 
defined aS the most nearly tangential orbits for which incident particles 
of a giiren energy are able to reach a given pOint on the spacecraft sur- 
face, having tangential velocity component In a given direction 
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to be directed inic & cavity, fisiieclaily severe ai^clng jircbleme ai-e khcwft to have 

occurred between electrcnlc cempcnente mounted inside e cavity at one ertd of the 
OSes spacecraft. 

Finally, we discuss some further implications of the multiple-root results 
shown in Tables 1 and 2 and Figure 5. Consider a situation involving two or more 
adjacent but isolated spacecraft surfaces which are made of the same material, 
and whose external conditions change with time, as in-the case of tlme-varying 
ambient distributions, or a Spacecraft rotation which carries these surfaces from 
sunlight Into shadow. Such a situation might Involve the continuous evolution of a 
single-root Into a multiple-root Situation, and the possibility Would then arise of 
a "bifurcation" phenomenon in which different surface elements followed different 
potential histories, with a Correspondingly large potential difference arising 

between them. Again, detailed numerical simulations are necessary in order to 
find out if stidh phenomena can actually occur. 


4, CONCLUSJOISS 

Upiief* bounds have been calcxilated for negative floating potentials tvhich may 
be acquired by electrically isolated Shaded Surfaces bn Synchronous Spacecraft. 
Effects of spacecraft shape and sheath potential profile haVe been Shown to have 
large Influences bn such potentials, inclusion of electron backs catterlng currents 
causes only a moderate reduction of these negative potentials in most cases. For 
isolated surfaces inside shaded cavities, negative floating potentials may exceed 
those on convex surfaces, m Some conditions, two possible floating potentials 
exist, leading to the possibility cf "bifurcation phenomena" in which adjaeeht isola- 
ted surfaces made of the same material may follow different charging histories. 
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